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How the extinction of extragalactic background light 
affects surface photometry of galaxies, groups and clusters 
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ABSTRACT 

The faint regions of galaxies, groups and clusters hold important clues about how 
these objects formed, and surface photometry at optical and near-infrared wavelengths 
represents a powerful tool for studying such structures. Here, we identify a hitherto 
unrecognized problem with this technique, related to how the night sky flux is typically 
measured and subtracted from astronomical images. While most of the sky flux comes 
from regions between the observer and the target object, a small fraction - the extra- 
galactic background light (EBL) - comes from behind. We argue that since this part of 
the sky flux can be subject to extinction by dust present in the galaxy/group/cluster 
studied, standard reduction procedures may lead to a systematic oversubtraction of 
the EBL. Even very small amounts of extinction can lead to spurious features in radial 
surface surface brightness profiles and colour maps of extended objects. We assess the 
likely impact of this effect on a number of topics in extragalactic astronomy where 
very deep surface photometry is currently attempted, including studies of stellar halos, 
starburst host galaxies, disc truncations and diffuse intragroup/intracluster light. We 
argue that EBL extinction may provide at least a partial explanation for the anoma- 
lously red colours reported for the halos of disc galaxies and the hosts of local starburst 
galaxies. EBL extinction effects also mimic truncations in discs with unusually high 
dust opacities, but are unlikely to be the cause of such features in general. Failure 
to account for EBL extinction can also give rise to a non-negligible underestimate 
of intragroup and intracluster light at the faintest surface brightness levels currently 
probed. Finally, we discuss how EBL extinction effects may be exploited to provide 
an independent constraint on the surface brightness of the EBL, using a combination 
of surface photometry and direct star counts. 

Key words: galaxies: photometry - diffuse radiation - dust, extinction - galaxies: 
spirals - galaxies: halos - galaxies: clusters: general 



1 INTRODUCTION 

The low surface brightness regions of galaxies, galaxy 
groups and clusters hold important clues about how such 
systems formed. The faint stellar halos of galaxies are 
predicted to contain signatures of p ast mergers in the 



hiera rc hical build-up of galaxies (e. g. Bullock fc Johnston! 



l2005l : lAbadi. Navarro fc Steinmet j |2006| ). and intra- 
group and intracluster light may help us understand 



NaDolitano et al.ll2004l:lMurante et al. 
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20061: Purcell. Bullock & Zenter 2007 


). The outermost re- 



gions of galactic discs may also allow us to constrain the 
prop ert ie s of star formation thresholds in the inter stellar 
medium (|Elmei2:reen fc Parravanolll994l : ISchavdl200j ). 

At the current time, there are two complementary tech- 



niques for studying the faint regions of extended objects: 
direct star counts and surface photometry. By resolving 
bright, individual stars, it is possible to trace the outskirts 
of galaxies to remarkable distances from their centres, cur- 
rently equivalent to surface brightness levels of /iy ~ 35 mag 
arcsec"'^ (e.g. [ibata ct al. 2007 ). However, this technique is 
only applicable for nearby systems. Moreover, it is currently 
only the brightest stars - with a very limited span of initial 
masses - that are probed, whereas more diffuse flux compo- 
nents produced by the myriad of faint stars below the detec- 
tion threshold and light emitted by the interstellar medium 
are not directly measured. Surface photometry offers a more 
concise census of the surface brightness, but is subject to a 
host of systematic problems which prevents regions as faint 
as those probed by star counts to be studied. 

The two chief problems which limit the depth of 
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optical /near-IR surfac e photometry are ins t rumental scat - 
tering llMichardl l2002l : ISirianni et all l2005l : Ide Jond |2008| ) 
and the challenge of subtracti ng the night sky flux with 
sufficient accuracy (e.g. Mclnick. Selman &: Quintanal 1 19991 : 
IZheng et all 1 19991 : iMihos et all 120051 ) . Because of these ob- 
stacles, optical surface photometry rarely probes regions 
faintward of the /x « 28 mag arcsec"^ isophote, although a 
few attempts have been made to push the limits into the /x ~ 
30-32 mag arcsec"^ range (e.g. iBarton fc Thompson 1 19971 : 



IZibetti. White fc Brinkmannll2004 : IZibetti et al.ll2005^ ■ De- 
spite these limitations, surface photometry is a technique 
that we will have to live with for a long time to come. 

Here, we identify a previously unrecognized problem 
with surface photometry at the faint limit, related to how 
the sky flux is typically estimated and subtracted from as- 
tronomical images. While most of the sky flux (in the form 
of airglow, zodiacal light, light from the stars and interstel- 
lar medium in the Milky Way) comes from regions between 
the observer and the object studied, a small fraction - the 
extragalactic background light (EBL) - comes from behind. 
The EBL at optical and infrared wavelengths is believed 
to be the product of direct and reprocessed starlight emit- 
ted over the entire star formation history of the Universe, 
and hence stems from objects at vastly different redshifts. 
Provided that the existing EBL measurements are correct, 
most of this light appears to be diffuse (i.e. unresolved with 
all existing instruments). Unlike the other components of 
the night sky flux, the EBL can be subject to extinction by 
dust present in the low surface brightness regions of the tar- 
get object, thereby invalidating an implicit assumption in all 
current surface photometry measurements, namely that the 
sky flux and the flux from the target object are unrelated. 

In Section 2, we explain how this effect is likely to man- 
ifest itself in the surface brightness proflles of extended, low- 
redshift objects, arguing that spurious features are expected 
to turn up at surface brightness levels similar to or fainter 
than that of the EBL itself. Estimates of the surface bright- 
ness of the EBL in units of mag arcsec"^ are presented in 
section 3 for Johnson- Cousins UBVRIJHK and Sloan Dig- 
ital Sky Survey (SDSS) ugriz filters. In sections 4 through 
7, we describe the relevance of EBL extinction for studies 
of disc truncations, halos of disc galaxies, host galaxies of 
local starbursts and intragroup/cluster light, respectively. 
In section 8, we outline how EBL extinction effects may be 
exploited to provide an independent measurement of the sur- 
face brightness of the EBL, using a combination of surface 
photometry and direct star counts. A number of remain- 
ing uncertainties in our analysis are discussed in section 9. 
Section 10 summarizes our findings. 



2 HOW EXTINCTION OF EXTRAGALACTIC 
BACKGROUND LIGHT AFFECTS SURFACE 
PHOTOMETRY 

The surface brightness of the night sky (/iv ~ 21.9 mag 
arcsec"^ in the y-band at the darkest telescope sites on 
Earth) is composed of several components: airglow, zodi- 
acal light, light from the Milky Way, and the EBL. The 
latter, which consists of the integrated contributions from 
astronomical light sources over a wide range of redshifts, 
contributes only a small fraction to the overall night sky 





Figure 1. Schematic illustration of how the surface brightness of 
the night sky is typically estimated. The sky level is measured in 
regions (boxes) well away from the target galaxy (central object) 
and other resolved sources in the frame. A surface or constant is 
then fitted to these measurements and subtracted from the frame. 
The problem with this procedure is that the sky level relevant for 
the target isophote (dashed ellipse) is likely to be slightly lower 
than in the regions where it was estimated. This happens because 
of EBL extinction by dust associated with the galaxy studied. 



flux (/^EBL,\^ ~ 25.4 arcsec"'^, as estimated in section 3, or 
about 4% of the total), but is the one component that is 
likely to originate largely from behind target objects in the 
low-redshift Universe. 

This turns out to be relevant when attempting to push 
surface photometry techniques to very faint surface bright- 
ness levels. To reach structures far below the surface bright- 
ness of the night sky, the subtraction of the sky flux needs 
to be as accurate as possible. As an example, surface pho- 
tometry measurments at ~ 28.0 mag arcsec"'^ with an 
uncertainty of ±1.0 mag require that the sky flux can be 
measured with an accuracy better than 99.9%, which is a 
daunting task at the current time. 

In Fig. [1] we illustrate how the sky flux is usually esti- 
mated from optical images. Consider the case of a galaxy, for 
which one would like to derive the surface brightness profile 
out to a radius marked by the dashed ellipse. One typically 
measures the sky flux in regions (here marked by boxes) well 
away from the target galaxy (at least outside the faintest 
isophote one aspires to reach) and other resolved sources 
in the frame. A surface or constant is then fitted to these 
measurements and subtracted from the frame. An implicit 
assumption in this procedure is that the sky level determined 
outside the target isophote (dashed ellipse) can somehow be 
interpolated into the regions where one attempts to measure 
the surface brightness of the galaxy studied. The problem 
is, that in the presence of dust within the target galaxy, the 
EBL component of the sky flux can be somewhat lower in- 
side the galaxy than outside. Hence, the sky flux measured 
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Figure 2. Schematic illustration of the spurious effects introduced in surface brightness profiles by failure to take EBL extinction into 
account, a) Typical sky subtraction procedures assume that the EBL (thin gray line), and hence the night sky (thick gray line), is 
featureless across the surface brightness profile of an extended object (black line), b) In reality, the surface brightness of the EBL may, 
due to extinction by dust present within the object, show a central depression across the face of the object, resulting in a central depression 
in the total surface brightness of the night sky. In relative terms, the depression in the total night sky is significantly weaker, due to 
dilution by other night sky components (airglow, zodiacal light, light from the Milky Way), which do not show any such depression, 
c) As a result of assuming the situation depicted in a), rather than b), current sky subtraction procedures are unlikely to recover the 
intrinsic surface brightness profile of the object (solid line), but instead a profile (dashed line) for which the surface brightness drops too 
fast at large radii. This will not have any significant effect on the total fiuxes and colours of bright objects, but may significantly affect 
studies of their low surface brightness regions. 



far away from the galaxy will systematically overestimate 
the sky flux at positions within the galaxy. 

The effect of this oversubtraction of sky flux on the 
resulting surface brightness profile of the galaxy is schemat- 
ically illustrated in Fig. [2] Fig. [2^ depicts the situation im- 
plicitly assumed when reducing astronomical images. The 
EBL (thin solid gray line) and consequently the total sky 
background (thick solid gray line) are incorrectly taken to be 
featureless across the surface brightness profile of the target 
object (solid black line). In Fig.[2jD, we depict what is likely 
to be a more realistic situation, in which the EBL displays 
a depression due to dust extinction in the direction towards 
the target object. As a result, the total sky background also 
shows a decrease in the direction towards the target. Fig. [2}; 
shows the resulting surface brightness profiles after sky sub- 
traction. The incorrect assumption of a flat sky background 
across the object gives rise to an oversubtraction of sky, and 
hence a surface brightness proflle which drops too fast at 
large distances from the centre (dashed black line) compared 
to the intrinsic profile (solid black line). Since the EBL con- 
stitutes such a small part of the overall sky flux, this problem 
will have negligible effects on the total and colours fluxes of 
bright galaxies, but may signiflcantly affect surface photom- 
etry measurements of low surface brightness regions. One 
may argue that such regions are unlikely to contain much 
dust, but we will demonstrate that even tiny amounts of ex- 
tinction can have significant effects at the faintest isophotes 
currently probed. 

The sky subtraction situation just described matches 
that typically encountered when using images built from a 
stack of short exposures taken with small telescope offsets 
inbetween (a technique sometimes referred to as dithering). 
In the near-IR, one often removes the sky from dithered 
images by subtracting subsequent exposures from each other 
(provided that the offset has been greater than the diameter 
of the target isophote). In the optical - and especially in the 
near-IR - one sometimes also choses to move the telescope to 
one or several adjacent sky fields inbetween exposures of the 
target object (sky chopping). Either the entire sky frame, or 
a surface fitted to it, may then be subtracted from the object 



frame. However, these alternative sky subtraction techniques 
are all based on the same basic assumption as the method 
first outUned, namely that the sky fiux outside the target 
object is representative for the sky flux inside. Therefore, 
all of these methods will result in a similar oversubtraction 
of the EBL. 

The effects of oversubtracting the EBL can be estimated 
using a set of simple equations. Standard reduction tech- 
niques (which neglect the effects of dust on the EBL) im- 
plicitly assume that the observed surface brightness I of the 
target (in units of Lq pc~^) at a certain wavelength A and 
position r is given by: 

-'obj, approx 

( A, r) = /obj ( A, r) + /airgiow ( A, r) + 

-^zodiacal 

/mw(A, r) + Jebl(A, r), (1) 

where Jobj is the intrinsic surface brightness of the target 
object and /airgiow, /zodiacal, Imw and /ebl are the contri- 
butions from airglow, zodiacal hght. Milky Way sources and 
EBL to the surface brightness of the night sky. However, 
since the EBL has passed through the dust layers of the 
target object, the observed surface brightness is more accu- 
rately given by: 

/obs(A,r) = Jobj (A, r) + 

-^airglow 

(A,r) -I- /zodiacal ( A, r) -I- 

/Mw(A,r) + (1 - /abs(A,r-))/EBL(A,r), (2) 

where /abs(A,r) corresponds to the fraction of the EBL flux 
that is absorbed by dust at wavelength A and position r. 
Under the incorrect assumption that equation ([ij holds, 
standard sky subtraction methods would result in a surface 
brightness /obj, approx(A,r) for the target object: 

/obj , approx — /obs 

(A,r) - /airglow (A, r) - 

/zodiacal (A, r) -/Mw(A,r-) -/EBL(A,r-), (3) 

whereas the true, intrinsic surface brightness of the object 
is given by: 

/obj(A,r) = /obs(A,r) — /airglow (A, r) - 

/zodiacal (A, r) - /Mw(A,r-) - (1 - /abs(A,r))/EBL(A,r). (4) 

Hence, because of the incorrect assumptions going into 
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equation ((TJ, standard sky subtraction methods introduce a 
relative error in the surface brightness of the target galaxy 
corresponding to: 



/obj(A,r) - /obj, 



^bj(A,r) 

/ab3(A,r)JEBL(A,r-) 

^bj(A,r) 



(5) 



For reasonable values on the dust absorption factor 
(/abs "C 1), this error will be small when the surface bright- 
ness of the target object is greater that that of the EBL 
(-'^obj(A,r) S> /EBL(A,r)). It may, however, become consid- 
erable when the surface brightness of the target becomes 
similar to, or lower, than the EBL (7obj(A,r) /EBL(A,r). 
Since the error introduced by the implicit assumptions going 
into equation ([T} results in a systematic oversubtraction of 
sky flux, this error is always positive. This means that the 
surface brightness of the target object will be systematically 
underestimated. In Fig. we demonstrate this effect on 
exponential surface brightness profiles 7(r) = /oexp(— r//i) 
(often used to describe the surface brightness profile of disc 
galaxies) in the case of several different spatially constant 
dust absorption factor /abs (A) (see Section 4 for a treatment 
of more realistic disc opacity profiles). The deviation of the 
observed surface brightness profile from the intrinsic one be- 
comes significant around Jobj ~ Iebl, and keeps getting in- 
creasingly serious at fainter surface brightness levels. Once 
/abB(A,r)7EBL(A,r) > Jobj, /obj, approx wiU bccomc negative, 
at which point the observer will supposedly assume that the 
signal has been lost in the sky noise. 

As seen in Fig. even an EBL extinction as low as 
A{V) = 0.01 mag (corresponding to /abs(V^) ~ 0.009) may 
give rise to detectable effects when pushing surface bright- 
ness measurements to their current limits (logjQ I{r) //ebl ^ 
—2, corresponding to fiv ^ 30.5 mag arcsec"^). 

The discrepancy between the true and inferred surface 
brightness has a complicated wavelength dependence, which 
may introduce spurious colour gradients across the face of 
extended objects. In the case of a colour index mA,i — mx.2 
(e.g. V — I, measured in magnitudes), the failure to account 
for EBL extinction will introduce a shift of 



A(mAi -mx^) = -2.51ogi 



Jobj(^l.'-) 



(6) 



We stress that this colour shift is different from nor- 
mal dust reddening, because of its dependence on the ratio 
/obj/^EBL- In principle, the EBL effect can shift the colour 
either redward or blueward depending on the extinction law, 
the spectrum of the target object and the spectrum of the 
EBL itself (see Sections 5 and 6 for examples). 

In Fig. Isjj, we use the values estimated in section 3 for 
the surface brightness of the diffuse EBL in filters V and 
I to demonstrate the spurious colour gradients introduced 
in V — I across an exponential disc when EBL extinction 
is neglected. At surface brightness levels comparable to or 
lower than that of the EBL, the observed V — I colour is 
driven away from the intrinsic value, eventually reaching 
values not expected for normal stellar populations. 



3 SURFACE BRIGHTNESS OF THE 
OPTIC AL/NEAR-IR EBL 

As illustrated in Fig. [T] it is customary to mask resolved ob- 
jects (foreground stars and background sources) before esti- 
mating the sky level to be subtracted from optical/near-IR 
images. As resolved background objects make up a fraction 
of the EBL, the part of the EBL that is relevant for EBL 
extinction effects is the diffuse component, i.e. the fraction 
of the EBL that remains after masking. We stress, however, 
that this component need not be diffuse in any absolute 
sense - the only requirement is that the sources responsible 
for the diffuse EBL remain unresolved (and unmasked) in 
the image analyzed. Deep, high-resolution images will there- 
fore contain a diffuse EBL smaller than that in shallow im- 
ages taken under poor seeing conditions. 

The surface brightness of the total EBL at optical/near- 
IR wavelengths remains very uncertain, mainly due 
to the problem of corre c ting f or zodiacal light (e.g . 
Dwek. Arendt fc KrennrichI l2005l: iMatsumoto et all l2005l : 



Mattilal l2006l : [Thompson et all l2007l : iBernsteinI |2007^ . hi 
what follo ws, we adopt the direct EBL measurements pre- 
sented by IBernsteinI (|2007| ) in the Hubble Space Telescope 
(HST) WFPC2 fihers F300W, F555W and F814W filter 
(similar to Johnson-Cousins (7, V a nd 7). In the near-IR, 
we adopt the J-band measuremen t of|W right (2001*) a nd the 
K and L band measurements of iGoriia n. Wright fc Charvl 
(|2000l ) . By linearly interpolating the EBL spectrum between 
the measured data points, we estimate the surface brightness 
of the total EBL in Johnson/ Cousins UBVRIJHK filters 
and SDSS ugriz filters. The result is presented in the sec- 
ond column of Table 1. While EBL measurements are often 
presented in units of ergs ~^ cm~^ sr~^ or similar, we 
here present all estimates in units of mag arcsec"^, which 
are more commonly used in the field of surface photometry. 
Vega magnitudes are used for UBVRIJHK and SDSS AB- 
magnitudes for ugriz. The errorbars presented simply reflect 
the errors quoted for the original EBL measurements, and 
do not take into account the errors i ntroduced by the inter - 
polation scheme. As discussed bv .Matsumoto et al.l l|2005h . 
the direct EBL measurements hint at a peak in the EBL 
spectrum around the J band (but see iMattilal |2006| and 
[Thompson et al.ll2007l for different views), and the unknown 
shape of this putative feature makes the a-band (interpola- 
tion from 7 to J) estimates for the EBL rather uncertain. 

While some of the EBL implied by these measurements 
can be attributed to resolved galaxies, m ost of it remains un- 
resolved even in the deep est 775r images (jPozzetti fc Madaul 
I2OOII : iTotani et al.ll200ll ). Indirect constraints based on the 
opacity of the Univers e at gamm a -ray wavelengths (e.g. 
Franccschini, Rodighicr o &c Vaccaril 20081: Mazin fc Raud 
[20071 : lAlbert et al.li2008i : iRazzague. Dermer fc Finkd |2008| ) 
tend to favour an EBL level closer to that attributable to 
the resolved source s, but such methods may have problems 
of th eir own (e.g. iKrennrich et al. I I2OO8I : iDe AngeUs et all 
l2009f ). To derive the surface brightness o f the resolved EBL, 
we adopt the estimates computed by iMadau fc Pozzettil 
(2000) by integrating HST galaxy counts in a set of 7 
broadband filters (somewhat similar to Johnson- Cousins 
UBRIJHK) from the Hubble Deep Fields. These galaxy 
counts are formally valid to faint AB magnitude limits vary- 
ing from 30.5 mag at 6700 A to 25.5 mag at 22000 A, but 
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Figure 3. The effects of spatially uniform EBL extinction (dust screen) on the V-band surface brightness profiles and colours of an 
exponential disc (central surface brightness /xq = 21 mag arcsec"^, scale length h = 2.5 kpc, intrinsic colour V — I = 1.2). a) The 
intrinsic V-band surface brightness profile (black solid line) compared to the profiles resulting from failure to account for EBL extinction 
(dashed lines). The gray horizontal line represents the EBL level. The different dashed curves correspond to (from left to right at 
logio I/Iebl = -2.5) extinction values of A{V) = 1.0, 0.1, 0.01 and 0.001 mag (/abs(^) ~ 0.6, 0.09, 0.009, 0.0009). The EBL extinction 
has pronounced effects on the profiles at surface brightness levels comparable to, or lower than that of the EBL itself. Even extinction 
values as low as A{V) = 0.01 mag may give rise to detectable effects when pushing surface brightness measurements to current limits 
(e.g. logiQ I /Iebl ~ —2, which corresponds to fiy = 30.5 mag arcsec"^ in the case of /^ebl, v = 25.5 mag arcsec" ^ as derived in 
Section 3). b) The spurious A(V — I) colour gradient introduced by EBL extinction in the case of A{V) = 0.1 mag and 0.01 mag as a 
function of intrinsic V-band surface brightness. In the case of A(V) = 0.1, significant colour offsets are introduced at log^Q I / Iebl < 0) 
(corresponding to > 25.5 mag arcsec" ^). For A{V) = 0.01 mag, the same effect is seen at log^Q I/Iebl < —1 (f^v > 28.0 mag 
arcsec" ■^). 



iMadau fc Pozzettil (|200Cll ) argue that the EBL contribution 
from resolved sources has converged in all the filters con- 
sidered. We caution, however, that alternative methods for 
estimating the resolved compone nt of the EBL give some- 
what different res u lts (compare iMadau fc Po zzctti 200Q : 
iTotani et al. 2001 : Bernstein et alT2002l ). By interpolating 
the iMadau fc Pozzettil (|2000l ) estimates, we arrive at the 
UBVRIJHK and ugriz estimates for the surface bright- 
ness of the resolved EBL presented in the third column of 
Table 1. 

To assess the surface brightness of the diffuse (unre- 
solved) EBL, we subtract the resolved EBL component from 
the total: 



L, diffusc(-^) — -^EBL, totzli^) — h 



EBL, resolved 



(A). 



(7) 



Here, we for simplicity assume that the all of the EBL that 
can currently be attributed to resolved galaxies has been 
resolved and masked in the images on which surface pho- 
tometry is attempted. While this is likely to hold only for 
deep images taken from space or under excellent seeing con- 
ditions from the ground, the corrections due to magnitude 
limit s different fro m those adopted above can be assessed 
from lPozzetti fc"M adau (2001,). Under the assumptions that 
the total EBL measurements cited above do not contain any 
major systematic errors, corrections of this type turn out 
to be modest in all filters considered. The resulting surface 
brightness of the diffuse EBL are presented in the fourth 
column of Table 1. We make no attempt to assess the error 
bars on these entries, since systematic effects are likely to 
dominate the uncertainties. 

fnspection of the entries in Table 1 reveals that the 
surface brightness of the total EBL lies in the range fi ~ 25- 



26 mag arcsec"'^ in the optical (UBVRI) and fj, « 23-25 
mag in the near-IR (JHK). This is considerably brighter 
than the isophotal levels typically probed through optical 
surface photometry, and around the limit of what ground- 
based near-IR surface photometry can currently reach. Since 
the resolved galaxy population seems to make up no more 
than « 1/10 (optical) to 1/3 (near-IR) of the total EBL, 
the surface brightness of the diffuse EBL is only 0.1-0.5 mag 
fainter than the total. 

Given that EBL extinction is likely to become relevant 
at surface brightness levels similar to or fainter than that of 
the diffuse EBL (see Fig. [3}, one may suspect this mecha- 
nism to impose strong limits on how deep surface photom- 
etry measurements can be pushed without the introduction 
of systematic errors. In Fig. 3] we show the predicted offset 

between the intrinsic and observed surface brightness, 
as a function of observed surface brightness /iobj, approx (i.e. 
the values inferred when EBL extinction is neglected) , in fil- 
ters UBVRIJHK for spatially constant V-band extinction 
values of A{V) = 1.0 mag (thick lines), 0.1 mag (medium 
lines) and 0.01 mag (thin lines), assuming a Milky Way ex- 
tinction curve. These offset predictions are independent of 
the surface brightness profile of the target object (as long as 
the extinction is considered spatially constant), and can be 
used to give a rough assessment of the likely impact of EBL 
extinction on any set of surface photometry data. The off- 
sets turn out to be fairly independent of the filters used. This 
"conspiracy" stems from the fact that the brighter EBL in 
the red part of the spectrum (when measured in Vega mag- 
nitudes) is compensated for by the lower extinction given by 
Milky Way extinction at those wavelengths. 

Systematic errors of ~ 1 mag are produced at 
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Table 1. The estimated surface brightness of the total, resolved and diffuse EBL in Johnson-Cousins UBVRIJHK and SDSS ugriz 
filters. The surface brightness of the diffuse component is obtained by subtracting the EBL attributable to currently resolved galaxies 
from the total EBL. UBVRIJHK magnitudes are in the Vega photometric system and ugriz magnitudes in the SDSS AB system. See 
main text for additional details. 



Filter Total EBL Resolved EBL Diffuse EBL 

(mag arcsec^^) (mag arcsec"^) (mag arcsec"'^) 



u 


26.0 


+0.9 
-0.5 


28.4 


+0.2 
-0.2 


26.1 


B 


26.1 


+0.S 
-0.4 


28.6 


+0.2 
-0.1 


26.3 


V 


25.4 


+0.7 
-0.4 


27.9 


+0.2 
-0.1 


25.5 


R 


25.0 


+0.8 
-0.4 


27.4 


+0.2 
-0.1 


25.1 


I 


24.6 


+0.9 
-0.5 


26.7 


+0.2 
-0.1 


24.7 


J 


24.3 


+0.5 
-0.4 


25.5 


+0.2 
-0.2 


24.7 


H 


23.9 


+0.3 
-0.3 


25.1 


+0.3 
-0.2 


24.4 


K 


23.1 


+0.3 
-0.3 


24.2 


+0.2 
-0.2 


23.5 


u 


26.7 


+0.9 

-0.5 


29.1 


+0.2 
-0.2 


26.8 


9 


25.8 


+0.7 
-0.4 


28.3 


+0.2 
-0.1 


25.9 


r 


25.2 


+0.7 
-0.4 


27.7 


+0.2 
-0.1 


25.4 


i 


25.0 


+0.8 
-0.5 


27.3 
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^ ~ 27 mag arcsec"^ for A{V) = 1.0 mag, at ^ f» 29 mag 
arcsec"^ for A{V) = 0.1 mag and at /i « 31.5 mag arcsec"^ 
for A(V) = 0.01 mag. Significant EBL effects are therefore 
expected in current optical surface photometry data (which 
can reach as deep as ~ 32 mag arcsec"'^), whereas only mod- 
est effects are expected in the near-IR bands (JHK), since 
current observations rarely probe beyond fi ~ 23-24 mag 
arcsec"^ at these wavelengths. 



4 IMPACT ON STUDIES OF DISC 
TRUNCATIONS 

Fig. 13] suggests that surface brightness, when studied 
through deep surface photometry, should fall faster than 
that for an exponential in the outer, faint regions of discs. 
Such truncated surface brightness profiles are in fact com- 
monly observed am ong disc galaxies (for a recent review, see 
iPohlen et aLll2004l ). and one might ask whether some frac- 
tion of these features may be due to EBL extinction effects 
rather than intrinsic disc truncations. The possibility that 
artificial disc truncations may arise from oversubtraction of 
sky ha s previously bee n discussed b vlBarteldrees fc Dettmail 
l|l994h Ide Griis et al] (I2001I ) and iNaravan fc Jod (|2003h . 
but the effect discussed here is markedly different, since 
the amount of oversubtraction associated with the EBL 
is a function of both intrinsic disc surface brightness and 
dust content. Observationally, the truncation radius of both 
high and low surface brightness discs appears at an aver- 
age J-band surface brightness of /x/ ~ 25.3 mag arcsec"^ 
jKregel fc van der Kruitl2004f ). Inspection of Table 1 reveals 
that this is fainter than the diffuse EBL level inferred by di- 
rect measurements {fij ~ 24.7 mag arcsec"^), and hence in 
the surface brightness regime where EBL extinction effects 
are expected to turn up. 

To assess the magnitude of the EBL effect in the out- 
skirts of discs, we adopt a model for radial disc opacity r(r) 
of the type: 



r(r) — r(0) exp ( 



(8) 



where r(0) is the central opacity and hr is an extinction scale 
length. This model has been demonstrated to give a rea- 
sonable fit to the observed opacity distribution, as derived 
both from occulting galaxy pairs and counts of background 
galaxies through foreground dis cs (e.g lHolwerda et al. I I2OO5I : 
iHolwerda. Keel fc Bohonll2007l ). In what follows, we adopt 
To — 0.2 and /ir = 2.9i?off (where R^s is the effective ra- 
dius of the dis c), as derived for 0.01 < z < 0.1 galaxies 
in the SDSS bv lHolwerda et al.1 |20o3). This opacity profile 
was derived by stacking data in the SDSS r and i filters. 
To derive the opacity at other wavelengths, we assume the 
original opacity profile to be valid at A ~ 6900 A and sim- 
ply re scale tq b y assuming either a Milky Way extinction 
curve (|Peilll992l ) or gray extinction. The latter option is sup- 
ported by absence of any obvious trend between the colours 
of background galaxies s een through foregroun d discs and 
the inferred disc opacity (|Holwerda et al.ll2005|) and can be 
attrib uted to a clumpy interstellar medium ( Keel fc White! 
l200lD . 

Using the relation J?efF ~ 1.68/i for an infinite exponen- 
tial disc, we get a ratio between the opacity scalelength and 
the scalelength of the luminous disc of h-r/h ~ 5. 

In Fig. O we demonstrate the effects of EBL extinction 
on a typic al disc galaxy with ixp.y = 21.0 mag arcsec"^, 
h = 4 kpc (|MacArthur. Courteau fc Holtzm an 2003) in fil- 
ters V and I, assuming intrinsic colours of V — I = 1.2, 
and no population gradient. As seen in Fig.[5K, the EBL ex- 
tinction gives rise to a smooth change in slope of the outer 
surface brightness profile (at around r J5 25-30 kpc or ~ 6-8 
scalelengths from the centre) . At this distance from the cen- 
tre, the extinction is predicted to be around 0.10-0.05 mag- 
nitudes (depending on the extinction law) in the l/-band 
and around 0.05 magnitudes in the /-band. This break re- 
sembles the truncations commonly seen in the outer parts of 
disc galaxies, but here occurs around fj.i ^ 27 mag arcsec"^ 
compared to the observed val ue of fii « 25.3 mag arcsec"^ 
(|Kregel fc van der Kruitll2004l ) . The choice of extinction law 
has but a minor impact on the appearance of this feature, as 
evident from the small difference between the thick dashed 
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Figure 4. The predicted offset Afi between the intrinsic and observed surface brightness, as a function of observed surface brightness /i 
in different filters for A(V) = 1.0 mag (thick lines), A{V) = 0.1 mag (medium lines) and A{V) = 0.01 mag (thin lines). A Milky Way 
extinction curve has been assumed, a) Filters U (solid), B (dashed), V (dash-doted) and R (dotted), b) Filters I (solid), J (dashed), H 
(dash-dotted) and K (dotted). 



line (Milky Way extinction) and thin dashed line (gray ex- 
tinction) in Fig. [5^. This implies that the EBL extinction 
effect cannot be the primary cause of the observed disc trun- 
cations. Since the rapid decrease in surface brightness caused 
by EBL extinction is predicted to set in at slightly different 
radii in the V and / bands, a strong, spurious colour gra- 
dient will develop in the outskirts of the disc (Fig. [S]d) . At 
r « 25-30 kpc, — / is suddenly driven from the intrinsic 
value of — / = 1.2 to — / > 1.6, which is not expected 
for normal stellar populations, even at high metallicities. 

While EBL extinction is unlikely to be the main cause 
of the observed disc truncations, the agreement between the 
observed and predicted disc breaks improves for dusty discs. 
This is demonstrated in Fig. (6] where the ratio between 
the opacity scalelength and the scalelength of the luminous 
disc has been increased to hr/h « 20 (in agreement with 
th e average value found for the smaller sample presented 
by iHolwerda et al] |2005| ) and the central opacity has been 
boosted to tq = 0.5. This places the spurious break in the I- 
band profile at around fj,i « 25 mag arcsec"^, in reasonable 
agreement with the break radius typically observed. 

This opacity model gives a central V-band EBL extinc- 
tion of A{V) ^ 0.7 (0.55) mag and A{V) ^ 0.55 (0.4) 
mag at 25 kpc from the centre for a Milky Way (gray) 
extinction curve. Assuming an intrinsic B — V colour of 
B — V = 0.6, this distance corresponds to r ~ 1.5 x R25 for 
the model disc, where 7?25 is the radius of the hb = 25 mag 
ar csec"^ isophote. Inspe ction of the opacity profiles derived 
in lHolwerda et al.1 Hooi) and lHolwerda et all (|2007| ) reveals 
that an opacity model of this type may well be consistent 
with some of the more opaque disc galaxies in their samples 
(although the error bars are typically very large for indi- 
vidual galaxies at these radii). This implies that while EBL 
extinction is unlikely to be the main reason for the observed 
disc truncations, some of the truncations seen in unusually 
dusty discs may possibly be attributed to this mechanism. 
While disc truncations often manifest themselves as a break 
between two exponential profiles with different scale lengths 
(in plots with /i versus R such as Fig. [5^ and Fig. |6l this 



would correspond to a straight line which suddenly changes 
slope at the break radius), EBL extinction effects predict a 
surface brightness profile which turns progressively steeper 
and eventually becomes vertical (see Fig. Fig. [S^ and 
Fig. [BJ. The distinction may, however, be difficult to make 
in the presence of large observational errorbars at the fainter 
isophotes, and the fact that surface photometry of discs 
rarely probe regions beyond = 28 mag arcsec"^. At this 
radius, the vertical part of the profile (where the surface 
brightness tends to zero due to oversubstraction of sky) has 
not yet been reached (see Fig. [5^ and Fig. |6} . Moreover, any 
deviation from the smooth opacity profile given by equation 
((HI will also distort the appearance of this feature. 

A strong colour gradient in the outskirts of discs, like 
the one leading to unphysically high values for V — I in 
Fig. [SJa, could in principle serve as a tell-tale signature of 
EBL extinction. However, inspection of equation ([6)l reveals 
that while similar colour gradients are expected for many 
combinations of EBL and galaxy spectra, they are by no 
means generic. As an example, consider the case where the 
spectrum of the EBL and the target galaxy have similar 
shape (/obj(Ai,r)//obj(A2,r-) « /ebl(Ai, r-)//EBL(A2, r)). In 
the presence of gray extinction (/abs(Ai,r) « /abs(A2,r)), 
one gets A(mAi —mx^) ~ 0, leading to no detectable colour 
gradient. Hence, while strange colour gradients in the outer 
regions of discs may serve as an indication that the EBL 
extinction effect is at work, the absence of such gradients 
does not necessarily prove that such effects are unimportant. 



5 IMPACT ON STUDIES OF THE HALOS OF 
DISC GALAXIES 

Attempts to probe the stellar halos of disc galaxies through 
surface photometry remains extremely challenging due to 
the very faint surface brightness levels at which such struc- 
tures are expected to become dominant (nv ^ 29 mag 
arcse c"^; lAbadi et al. l l2006l : lGauthier. Dubinski fc Widrowl 
I2OO6I ). Detections of such halos around edge-on disc galaxies 
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Figure 5. The predicted eflfects of EBL extinction on an exponential disc witli central surface brightness fiv.o = 21 mag arcsec"^, scale 
length h = 4 kpc, an intrinsic colour of V — I = 1.2 at all radii and opacity properties as specified in the main text, a) The intrinsic 
/-band surface brightness profile (thick solid) compared to the observed profiles predicted for a Milky Way extinction curve (thick dashed) 
and gray extinction (thin dashed). The horizontal, gray solid lines indicate the EBL level in this filter. While the EBL extinction effect 
mimics a soft disc truncation, the c hange in profile slope occurs a t much fainter isophotes {fij « 27 mag arcsec"'^) than for typical disc 
galaxies (fij 25.3 mag arcsec~^: iKregel &: van der KruitI |2004| ) . b) The intrinsic V — I colour profile (thick solid) compared to the 
predicted colour profile in the case of a Milky Way extinction curve (thick dashed) and gray extinction (thin dashed). Due to a more 
rapid decrease in V-band fiux at 25-30 kpc from the centre compared to the 7-band, V — / is driven to extremely high values at these 
radii. 




20 
R (kpc) 

Figure 6. Same as Fig. [5^, but for a disc with r = 0.5 and hr/h = 
20 instead of instead of r = 0.2 and hr/h = 5. The boosted 
opacity of this disc shifts the drop in the surface brightness profile 
to brighter isophotes (fij 25 mag arcsec"^ instead of fij ^ 27 
mag arcsec"^ for the more transparent disc adopted in Fig. lEK) 
and smaller radii. 



have been claimed (e.E 


t. Sackett et al.''l994': 'Lcaueux et al. 


199e; 


Rudv et all 1 1997 


: Ijames & Casalii]1998i; .Zibetti et al. 


2004; 


Zibetti & FerEuson 


12004: ICaldwell & BerEvalll 2006h. 



but many of these remain controversial due to pr oblems in 



subtr acting the sky with sufficient accuarcy iZheng et al 



199£ ) and correcting for instrumental scattering (|de Jone 
20081 ). A common feature of the reported halos is a very red 
spectral energy distribution, manifesting itself as excess flux 
with respect to expectations for normal stellar populations 
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Figure 7. The intrinsic g-band surface brightness profile (thick 
solid) compared to the observed profiles (thick dashed) predicted 
for a spatially constant extinction of A{g) = 0.1, 0.05 or 0.01 mag 
(from left to right at fig = 32 mag arcsec"^). The horizontal, gray 
solid line indicates the EBL level in this filter. 



in the wavelength range from i 7600 A) to K (« 22000 
A). If correct, this could be indicative of a st ellar halo with 
a very bottom-heavy in i tial mass function jLegueux et al.l 
1 19961 : iRudv et al.l Il997l : IZackrisson et al.l |2006| ). with im- 
portant implications for the missing-baryon problem (e.g. 
IZackr isson fc FIvnn! l200l ). 

Figs. |3] and [5] have already indicated that EBL extinc- 
tion may produce artiflcial reddening of stellar populations 
at faint surface brightness levels. Hence, it becomes relevant 
to explore whether this could be the reason for the anoma- 
lous colours reported for the halos of disc galaxies. We here 
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focus on the halo measureme nts made in the SPSS gri filter s 
(used bv lZibetti eraLll2004l and lCaldwell fc Bergvalill2006h . 
In Fig. [T] we demonstrate the predicted effects of EBL ex- 
tinction on the (/-band data of a halo following the Sersic 
profile: 



I = 7eff exp (-6n [(r/reff)^^" ~ ) > 
with bn ^ 2n~ 0.324 (jCiottilll99ll ). 



(9) 



n — 7, rcs ~ 10 kpc. 
This intrinsic profile has been adopted for illustrative pur- 
poses only, with the values f or n and reg chosen t o allow the 
surface photometry data of IZibetti et al.l l)2004h to be ap- 
proximately reproduced, after EBL extinction effects have 
been introduced. However, our results concerning EBL ex- 
tinction as the source of abnormal halo colours do not de- 
pend on the details of the intrinsic profile. 

The thick solid line in Fig.[7]corresponds to the intrinsic 
surface brightness profile whereas the dashed lines indicate 
the profiles produced by a spatially constant extinction of 
A{g) = 0.1, 0.05 or 0.01 mag (from left to right at fig = 32 
mag arcsec"^). With this amount of dust extinction, the 
observed surface brightness profile becomes severely affected 
a couple of magnitudes faintward of the diffuse EBL level 
(gray horizontal line), just like in the case of disc-like profiles 
(see Figs.|3]&[5l). 

Fig-EK indicates how the r — i vs. g ~r colours of an old, 
metal -poor (Z = 0.001) stellar halo population modelled 
using iMarigo et al] (|2008l ) isochrones with a Salpeter initial 
mass function would be affected in the presence of A{g) = 
0.1 mag extinction and a Milky Way extinction curve. The 
halo colours (cross indicating la error bars) reported by 
IZibetti et al. are far redder in r — i than the stellar 

population model at all ages, but EBL extinction effects 
may introduce radial colour gradients (dashed lines) which 
substantially improve the agreement between the model and 
observations. 

Fig. reveals a number of interesting features of the 
EBL extinction effect. Whereas normal dust reddening of a 
stellar population would shift the colours along a straight 
vector in a colour-colour diagram like this, EBL extinction 
shifts the colours along curved trajectories. Moreover, the 
length and direction of the colour shift depends on the in- 
trinsic colours of the affected stellar population. This is illus- 
trated by the three different offset vectors, originating from 
slightly different positions along the evolutionary track of 
the model stellar population. The reason for this behaviour 
is the dependence of the colour shift on the ratio between 
the galaxy and EBL fluxes in equation ((Gjl - if the colours of 
the galaxy change, so will the colour shift. Similar effects can 
also be produced by altering the colours of the EBL itself. 
The triangles along each offset track (dashed line) indicate 
the g-band surface brightness at which a particular observed 
colour has been reached. 

While EBL extinction clearly can reproduce the ob- 
served halo colours, this happens at a surface brightness 
level (fig « 30 mag arcsec"^) significantly fainter than that 



at which the plotted colours have been measured (/ig 



28 



the predictions of EBL extinction effects, one would have 
to boost the amount of extinction in the halo or alter the 
spectral energy distribution of the EBL. Adopting a gray ex- 
tinction curve also alters the colour shift, but typically not 
in the desired direction. Our tests indicate that, to repro- 



duce the halo colours at fig « 28 mag arcsec"^ while keeping 
the surface brightness of the EBL intact would require an 
EBL ext inction of A(g ) = 1.0 mag in the halo. Given that 
the iZibetti et al.l l|2004) halo colours are derived in a region 
located at a substantial distance away from the edge-on disc 
(a projected distance of twice the disc scale length), this 
amount of extinction seems unrealistically high. 

However, the halo colours can in principle be repro- 
duced with an EBL extinction of A{g) — 0.2 mag, if the 
total surface brightness of the EBL in the r and i bands 
is assumed to lie at the bright end of what is allowed by 
the error bars (this boosts the total EBL to fir ~ 24.8 mag 
arcsec"'^ and fii — 24.5 mag arcsec"^ and the diffuse EBL 
to fir = 24.9 mag arcsec"^ and fii — 24.6 mag arcsec"^). 
This is demonstrated in Fig. [8]3. Assuming the dust to be 
smoothly and uniformly distributed, and the halo stars to be 
isotropically distributed around the halo centre, one would 
expect the extinction optical depth experienced by the halo 
stars to be approximately half of that suffered by the EBL, 
meaning that A{g) = 0.2 mag for the EBL corresponds to 
A{g) = 0.1 for the halo stars, or a normal dust reddening 
of only E{g — r) « 0.03 mag and E{r — i) « 0.02 mag (too 
small to be plotted in Fig. [S| . 

As far as we know, there are no observational con- 
straints that reliably rule out this amount of dust extinction 
in the halo. Dust is well kno wn to exist a few kpc above 
the plane of edge-on discs (e.g. lAlton et al.ll200d l. but con- 
straints on the the opacity as a function of extraplanar dis- 
tance are scarce. By studying the reddening of background 
galaxies, Zaritsky (1994) claimed a tentative detection of 
A{B) f» 0.01 mag of extinction at 60 kpc into the halo 
of two nearby disc galaxies. More recently, I Menard et al.l 
(|2009h studied the reddening of quasars projected close to 
galaxies in the SDSS and derived an opacity profile for galac- 
tic halos from projected distances of 20 kpc out to several 
Mpc. At the innermost data point (« 20 kpc), the measured 
opacity was no more th an A{V) ~ 0. 3 ma g. Considering 
that the halo colours of IZibetti et al.l (|2004l ) were derived 
~ 8 kpc into the halo (based on the median scale length for 
galaxies in their sample), the extinction at that projected 
distance may of course be considerably higher. One would, 
however, naively expect the opacity of the halo to be lower 
than in a disc, and given the model described in section 4, 
the typical disc opacity at this surface brightness level would 
be A{g) < 0.1 mag, i.e somewhat lower than required to at- 
tribute the red excess of the halo solely to EBL extinction. 

Our treatment admittedly neglects the effects of instru- 
mental scattering on the observed surface brightness profile 
of the halo. Due to this effect, light from the central disc will 
be scattered outward and contaminate the halo region. The 
wavelength dependence of the instrumental point spread 
function (PSF) may by its elf introduce spurious colour gra- 
dients an d ide Jong (|2008|') ar gues that the halo signal re- 
ported by IZibetti et al.l ~ 20041 ) may be severely affected by 
this. If this is the case, then the light observed at any given 
radius in the halo will consist of at least two components: 
scattered light from the inner regions of the galaxy, where 
the dust content is likely to be higher but where the high 
galaxy-to-EBL fiux ratio ensures essentially no EBL extinc- 
tion effects, and very faint (fainter than observed) halo light 
which has suffered more EBL extinction effects than sug- 
gested by the observed surface brightness. Without a de- 
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Figure 8. a) The effect of EBL extinction in a diagram of r — i vs. g — r. The thick solid line indicates the evolution of a metal-poor 
(Z = 0.001) stellar population (thick solid) with a Salpete r initia l mass function and an exponentially decaying star formation rate 
(SFR oc exp—t/r with t = 1 Gyr), based on iMarigo et all (2008') isochrones. White ma rkers represen t ages of 1 Gyr (circle) and 14 
Gyr (star) . The cross represents the la error bars on the colours of the halo detected by IZibetti et al.l l|2004l ') around stacked edge-on 
discs in the SDSS (at fig Ki 28 mag arcsec"^). Dashed lines indicate how a spatially constant EBL extinction of A{g) = 0.1 mag would 
shift the observed colours away from the model predictions under the assumption of a Milky Way extinction curve. The differently-sized 
triangles along these tracks correspond to surface brightness levels of fig = 27, 28, 29 and 30 mag arcsec"-^ (largest to smallest triangles). 
EBL extinction can indeed produce colours close to those observed, but only at relatively faint surface brightness levels [fig ft! 30 mag 
arcsec"^). b) Same as a) but for Ag = 0.2 mag of EBL extinction and a diffuse EBL boosted to fir = 24.9 mag arcsec"^ and fii = 24.6 
mag arcsec"^ (at the bright end of what is allowed by the error bars on the total EBL in Table 1). In this case, t he obse rved halo can 



be approximately reproduced at the surface brightness level fig ft 28 mag arcsec 
The age of the halo population affected by EBL extinction is here set to 9 Gyr. 



at which they were reported bv iZibetti et al.i (I2003) . 



tailed model of the eflective PSF, it becomes very difficult 
to reliably predict how this will alter the eflects of EBL ex- 
tinction. 



One could argue that, if EBL extinction is the cause 
of the anomalously red halos of edge-on disc galaxies, then 
strange colours should also be seen at a similar surface 
brightness level in the outer parts of discs themselves. In 
fact, one would expect the dust extinction to be higher in 
the disc than in the halo, thereby making the spurious colour 
gradients along the disc even more pronounced. While it 
would be worthwhile to search for this effect, one should 
bear in mind that the disc is likely to have intrinsic colours 
that are different from the halo (thereby altering the colour 
offset vector, see Fig. |8K). The structure of the interstellar 
medium may also be different in the disc and halo, with 
implications for the grayness of the effective attenuation. 



As seen in Figs.|8^ and b, EBL extinction predicts that 
the colours should become more and more extreme as one 
approaches fainter and fainter isophotes. The detection of a 
radially steady halo colour would therefore speak against the 
idea of EBL extinction as the main cause of the red excess, 
unless one assumes a radial evolution in the dust reddening 
properties which happens to cancel the colour gradient. Of 
course, the intrinsic halo colour may not be radially con- 
stant, since stellar halos are likely to contain age and metal- 
Ucity gradients at some level. 



6 IMPACT ON STUDIES OF THE HOSTS OF 
LOW-REDSHIFT STARBURST GALAXIES 

Blue compact galaxies (BCGs) are metal-poor, low-mass 
syst ems currently undergoi ng intense star formation activity 
(e.g lKunth fc OstlinI [2OO0I ') . While the progenitors of these 
objects remain unidentified, deep surface photometry of the 
regions outside their bright star-forming centres have pro- 
vided important clues on the properties of the underly- 
ing s tellar population (e.g. iNoeske et al.l I2OO3I : ICaon etaP 
'2005!). Combined optical an d near-IR surfa ce photometry 
( BergvaU fc Ostlin 2002; Bcr gvall et al.ll200d ) have for some 
of these systems indicated colours that become progres- 
sively redder as one traces the surface brightness profile 
outwards, eventually reaching values that are difficult to 
reconcile with the type of stellar population one expects 
BCG host galaxies to have. This may be a signature of 
a host gal axy with a metallicity higher than the centra l 
starburst (jBergvall fc OstlinI liooj : IZackrisson et al.ll2009l ). 
or po ssibly a halo with a b ottom-heavy initial mass func- 
tion jZackrisson et all |2006^) similar to that advocated to 
explain the colours reported for the halos of edge-on disc 
galaxies. However, since these extreme colours turn up at 
surface brightnesses (fiy ~ 26-27 mag arcsec"^; fix ~ 22- 
23 mag arcsec"^) comparable to or below the EBL level, we 
here investigate whether EBL extincition could possibly be 
responsible. 

In Fig. EK, we compare the V — K vs. B — V colours of 
the f aint outskirts of BCGs (as presented bv lBergvall et all 
I2OO5I ) to the predicted evolution of a Z = 0.001 (close to the 
typical metallicity of ft; 10% solar measured from emission- 
line ratios towards the central regions of BCGs) stellar popu- 
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Figure 9. a) Observed V — K and B — V colours of the faint outskirts of BCGs (crosses indicating la error bars), compared to 
model predictions for a stellar population with metallicity Z = 0.001 and a Salpeter initial mass function, based on iMarigo et all feOOSi ') 
isochrones. White markers represent ages of 1 Gyr (circle) and 14 Gyr (star). The data points are clearly displaced to the right of the model 
tracks, indicating a X-band excess in the spectra of these objects. An exponentially declining star formation rate (SFR(t) oc exp—t/r) 
with T = 1 Gyr has been assumed, but using a different star formation history does not remove the apparent red excess, b) Same, but 
with error bars on the observed colours omitted to avoid cluttering. The dashed lines indicate how EBL extinction would shift the colours 
away from the model track in the case of A{V) = 0.5 mag and a Milky Way extinction curve. The differently-sized triangles along these 
tracks correspond to surface brightness levels of fiy = 26, 27 and 28 mag arcsec"'^ (largest to smallest triangles). For the chocen opacity, 
EBL extinction can in principle explain the observed colours at /ly ?s 26—27. The arrow in the lower part of the figure indicates the 
normal dust reddening vector resulting from this scenario. 



lation, based on lMarigo et al.l ()2008[ ) isochrones. A Salpeter 
initial mass function and an exponentially decaying star for- 
mation rate (SFR(t) oc exp(— t/r) with t — 1 Gyr) has 
been assumed. Many of the data points are 0.5-1.0 mag 
redder in V — if than the model tracks, indicating either a 
higher host metallicity, non-negligible dust reddening or an 
extreme initial mass function (for a recent discussion, see 
IZackrisson et al.ll2009l ). In Fig. |9Jd, we demonstrate how an 
EBL extinction of A{V) — 0.5 mag would shift the colours 
away from the Z — 0.001 track, under the assumption of a 
Milky Way extinction law. In this diagram, the error bars on 
the observed colours have been omitted to avoid cluttering. 
While the direction and magnitude of the shift is sensitive 
to the intrinsic colours of the stellar population (analogously 
to the situation in Fig.[S^), most of the observed colours can 
in principle be explained as due to EBL extinction effects. 
The one exception is the object a.t V — K ~ A: (Haro 11), 
but new data on this galaxy indicates a problem with the 
original measurement, and th at the redness of its outskirts 
may have been overestimated l|Micheva et al.ll2009l ). 

Assuming that the EBL experiences an opacity twice 
that of the stars in the host, an EBL extinction of AiV) = 
0.5 mag would correspond to normal dust reddening values 
oiE{B-V) « 0.08 mag and E{V-K) 0.22 mag (this red- 
dening vector has been indicated with an arrow in the lower 
part of Fig pj). This is withiii the lim it of E{B -V) < 0.1 
mag set by Bergvall fc OstlinI (|2002l ) through observations 
of the Hq/H/3 emission- line ratio towards the central star- 
bursts of some of these objects. While there are no direct 
constraints on the dust extinction in BCGs at large pro- 
jected distances from their centres, one would expect the 
extinction to be lower in the outskirts of these objects. At 
a surface brightness similar to that of the region where the 
colours are measured {^v ~ 27 mag arcsec"^), the disc opac- 



ity model adopted in section 4 would for instance give an 
EBL extinction of only AiV) = 0.1 mag. On the other hand, 
these are starburst galaxies, and the dust processing cycle 
may not behave the same way as in quiescent discs. Advo- 
cating E{B — V)^ 0.08 mag for the BCG host galaxy would 
certainly not seem any more controversial than assuming a 
metallicity higher than observed, or a bottom-heavy initial 
mass function. Hence, we conclude that EBL extinction re- 
mains a viable explanation for the red excess seen in the 
outskirts of BCGs. 



IMPACT ON STUDIES OF INTRAGROUP 
AND INTRACLUSTER LIGHT 



Optical surface photometry of galaxy groups 



White et al.l |2003|; lOa Rocha fc Mendes de Oliveiral 



Da Rocha et al.ll2008l) and clus t ers (e . g. iMihos et al] 



20051: 



2005 



2005 



Gonzales, Zabludoff & Zaritskv| l2005l : IZibetti et al 
Soigar, Graham & Joricn 2007|) have revealed a faint in- 
tergalactic flux component, stemming from stars liberated 
from tidally shredd ed galaxies in these dense environments 
l|Purcell et aLll2007h . Since current measurements reach as 
deep as fiv ~ 28.5 for individual systems (|Mihos e t al."2005|) 
and U p « 32 mag arcsec"^ for stacked clusters ([Zibctti ct cd] 
2005!), EBL extinction effects may well be an issue. As 
demonstrated in previous sections, the importance of such 
effects strongly depend on the amount of dust present at 
these faint surface brightness levels. Significant dust redden- 
ing of background galaxies have been reported for individual 
systems - most notably the M81 group (corresponding to 
A(V) K. 0.4 m ag if a Milky way extinction curve is assumed; 
IXilouris et al.l [2OO6), where the starburst galaxy M82 is 
ejecting dust into the intragroup medium. However, several 
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Figure 10. The offset A/i between the intrinsic and observed 
surface brightness, as a function of observed surface brightness 
IX, in filters g (solid hnes), r (dashed lines) and i (dash-dotted 
lines) predicted by EBL extinction effects over the range of surface 
brightnesses relevant for current measurements of intragroup and 
intracluster light. Line thicknesses indicate EBL extinction values 
of A(V) = 0.05 mag (thick lines), 0.01 mag (medium lines) and 
0.005 mag (thin lines), assuming a Milky Way extinction curve. 



investigations aiming to measure the reddening of galaxies 
and quasars behind large samples of clusters have reported 
no or very little average reddening ( A(V) ^ .06 mag, as- 
suming Milky Way extinction ; Chelou ch e. Koest cr & Bowen 
l2007l : lMuller et al.ll2008l : IBovv. Hogg fc Moustakasll2008l ). A 
highly irregular dust distribution may in principle give rise 
to sizable extinction even if there is lit tle reddening (as ap- 
pears to be the case in disc galaxies, see lKeel fc Whitg|200ll : 
[Holwcrda et al. 2005), but without detailed constraints on 
the grayness of the attenuation, we are forced to assume a 
Milky Way extinction law in the following. 

In Fig. 1101 we show the offset Afi between the intrin- 
sic and observed surface brightness, as a function of ob- 
served surface brightness fj, (i.e. the surface brightness in- 
ferred when EBL extinction is neglected), in filters g (solid 
lines), r (dashed lines) and i (dash-dotted Hnes) predicted 
for an EBL extinction of A{V) — 0.05 mag (thick lines), 0.01 
mag (medium lines) and 0.005 mag (thin lines). The differ- 
ence between Fig. [TO] and Fig. His that Fig. [TO] uses SDSS 
filters and dust opacities more relevant for the intracluster 
medium. 

In the case of A{V) = 0.05 mag (thick lines), the ob- 
served g-band profile is ~ 1.3 mag fainter than the intrinsic 
one at an observed fig ~ 30 mag arcsec"^ (approximately 
corresponding to the faint limit of the total o bserved light 
profi le derived from stacked cluster images bv lZibetti et al.l 
I2OO5I ). Even at fig « 28 mag arcsec (approximately the 
current limit of detections of intracluster light in individual 
clusters), the offset is « 0.3 mag. The effect becomes smaller 
for smaller A{V), and becomes irrelevant (< 0.05 mag) for 
A{V) — 0.005 mag at fig « 28 mag arcsec"^. 



8 A NEW ROUTE TO CONSTRAINING THE 
SURFACE BRIGHTNESS OF THE 
EXTRAGALACTIC BACKGROUND 

While EBL extinction appears to set a limit to how deep 
surface photometry can be pushed without introducing se- 
rious systematic errors, this effect need not only be a source 
of nuisance. As described in the introduction, there are two 
complementing methods for studying regions of low surface 
brightness: surface photometry and direct star counts. The 
radius out to which star counts can be trusted is primar- 
ily set by the level of contamination from foreground stars 
and background galaxies, rather than the brightness of the 
night sky. Hence, EBL extinction effects should not be an 
issue for surface brightness profiles based on this technique. 
This allows for an independent way of measuring (or set- 
ting an upper limit on) the surface brightness of the EBL, 
which could help resolve the controversy regarding current 
measurements (as described in section 3.). 

If surface photometry reveals a significant drop in sur- 
face brightness (as in Fig. |3^, Fig. [5^, Fig. |6] and Fig. [T]), 
whereas star counts indicate no such feature, this would be 
a tell-tale sign of EBL extinction effects. By looking for such 
discrepancies in systems with known opacity profiles, the 
surface brightness of the EBL in the filter used can directly 
be estimated from equation ((5)1. If, on the other hand, no 
such discrepancies between the two methods are detected, 
this can be used to set an upper limit on the EBL. 

The best targets for an endeavour of this type would 
seem to be disc galaxies, both because these have rela- 
tively high opacities (e.g. iHolwerda et al.l |2005| ) and also 
since average opacity profiles have already been derived (e.g 
iHolwerda et al.1 12007| ). Under the assumption of the disc 
galaxy model used in Fig. [S] the absence of any discrep- 
ancy between the surface brightness profiles derived from 
surface photometry and star counts (at the ^0.1 mag level) 
out to fiv ~ 28 mag arcsec"'^ would indicate that the dif- 
fuse V-band surface brightness of the diffuse EBL is a fac- 
tor of 5 lower (i.e. 1.75 mag fainter) than indicated by the 
direct measurements (Table 1). By combining surface pho- 
tometry, star counts and opacity measurements (through 
counts of background galaxies), a test of this type can in 
principle be implemented. The weakest link is likely to be 
the opacity measurements, since current techniques tend 
to yield ra ther large error bars on the extinction at faint 
isophotes (|Holwerda et al.l [200^. Combined surface pho- 
tometry and star counts are already available f or a small 
number of disc galaxies, f or instance NGC 4244 l|Frv et al.l 
l2008l : Ide Jong et al.1 [2007I ). but the edge-on orientation of 
this particular object makes it far from optimal. An alter- 
native strategy could be to analyze a sizable number of discs 
without indivi dual opacity mea s urem ents, since the extinc- 
tion profile of IHolwerda et al.l (|2007l ) is likely to hold on 
average. 



9 DISCUSSION 

9.1 The surface brightness of the diffuse EBL 

The quantitative EBL extinction effects discussed in previ- 
ous sections are based on the surface brightness of the diffuse 
(spatially unresolved) EBL, as estimated in section 3. Since 
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the diffuse EBL is derived by subtracting the resolved EBL 
from the total, we caution that the surface brightness of this 
component may change dramatically, if either the total EBL 
turns out to be overestimated (as implied by estimates of 
the opacity of the Universe at gamma-ray w avelengths; e.g. 
Franceschini. Rodighiero fc Vaccaril l2008l: iMazin fc Raug 



20071 : [Albert et al.N2008l : iRazzague. Dermer fc Finkd l2008l ) 

or, alternatively, if the EBL that can be attributed to re- 
solved galaxies has been underestimated. 

As shown in section 3, the current difference between 
the total and diffuse EBL is only 0.1-0.5 mag. This means 
that the contribution from resolved sources is minor and that 
EBL extinction effects are fairly insensitive to the masking 
of background galaxies. This may no longer be the case if 
the total EBL level is significantly lower than implied by 
current direct measurements. In principle, the diffuse com- 
ponent may tend towards zero if future measurements are 
able to completely reconcile the total EBL with resolved 
galaxies. We point out, however, that a too rigorous masking 
of faint background sources may give rise to unwanted ef- 
fects on its own. Due to contrast issues, background sources 
will be more difficult to identify when superposed on faint 
isophotes of the target object than when located at greater 
projected distances in the sky. Hence, it may be impossible 
to reliably mask out interlopers projected against the tar- 
get object as faint as those that can be masked in sky fields. 
The EBL level estimated after masking resolved sources well 
away from the target galaxy will thereby be too low, lead- 
ing to an undersubtraction of sky. This effect may to some 
extent be canceled by extinction within the target object, 
which makes background objects appear dimmer, but to es- 
timate the net result on colour maps and surface brightness 
profiles would require very careful modelling. 



9.2 Spatial fluctuations in the EBL 

Significant spatial fluctuations in the EBL have been 
reported on scales of ^ 1" and upw ards a t 1-2 
um (e.g. iKashlinskv et al.l |2002|: [O dcnwal d et al.l l2003l : 
iMatsumoto et aLlbOOsT iThomtason et al.ii2007i l. This means 
that the effects of EBL extinction are likely to depend on 
the angular size of the objects studied, since the EBL con- 
tribution to the overall sky flux may vary from target to 
target or across the face of a single object. This effect is 
difficult to quantify without knowledge of the detailed an- 
gular fluctuation spectrum as a function of wavelength, but 
it should be kept in mind that this may introduce significant 
uncertainties in the predictions for EBL extinction effects for 
individual objects. 



9.3 Impact of EBL extinction on studies of 
high-redshift galaxies 

As discussed in section 3, the discrepancy between direct 
EBL measurements and the estimated contribution from re- 
solved galaxies indicates that most of the EBL is diffuse (or 
currently unresolved). Current models moreover have great 
difficulties in reproducing the observed properties of the 
opti cal/near-IR EBL without resorting to exo tic scenarios 
f e.g. IMatsumoto et al.|[2005l : lFardal et al.ll2007h . Hence, the 
origin of the EBL and its redshift evolution remains poorly 



constrained. This makes it very difficult to assess the effect 
of EBL extinction on surface photometry of objects beyond 
the local Universe, since part of the EBL may then con- 
stitute a foreground, rather than a background, and would 
not give rise to effects such as those discussed in previous 
sections. 



9.4 Impact of EBL extinction on studies of 
point-like objects 

So far, we have discussed the impact of EBL extinction 
on surface photometry of extended objects, but faint point 
sources could in principle also be subject to similar effects. 
If the part of the EBL that current telescopes are unable 
to resolve is due to faint, extended wings of distinct astro- 
nomical objects or to some unknown mechanism in the in- 
tergalactic medium (i.e. a truly diffuse EBL), point sources 
with dust structure on the scale of the seeing disc could be 
affected in a way similar to that of extended objects. Once 
the surface brightness across the seeing disc becomes similar 
or lower than that of the EBL, an overdensity of dust corre- 
lated with this seeing disc (dust either in the direct vicinity 
of the point source or in the foreground) would make esti- 
mates of the sky flux (based on measurements outside the 
seeing disc) too high. This results in an oversubtraction of 
sky, in direct analogy with the effect discussed in previous 
sections. At such small angular scales, the issue of spatial 
variations in the EBL of course becomes severe, making it 
very difficult to assess the likely impact of EBL extinction 
on realistic sources. 



9.5 Impact of EBL extinction on studies of Milky 
Way objects 

This paper has focused on the impact of EBL extinction on 
surface photometry of extragalactic objects, but this mech- 
anism should also affect surface photometry measurements 
on more nearby objects, such as star clusters and nebulae 
within the Milky Way. The resulting oversubtraction of sky 
may in this case be even more severe, since the EBL does 
not necessarily represent the only sky component coming 
from behind the target objects - faint Milky Way stars or 
diffuse light from the background interstellar medium may 
also contribute. 



10 SUMMARY 

Our conclusions can be summarized as follows: 

(i) We have pointed out a previously unrecognized prob- 
lem with surface photometry at the faint limit, related to 
how the sky flux is typically estimated and subtracted from 
astronomical images. Since the EBL part of the sky flux is 
likely to originate from behind extended objects in the low- 
redshift Universe, the overall sky flux across such targets 
may be diminshed due to extinction by dust. This causes 
a systematic oversubtraction of sky, which will affect sur- 
face brightnesses similar to or fainter than that of the dif- 
fuse (unresolved) component of the EBL. Without detailed 
knowledge of the opacities of the objects studied, this effect 
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sets a limit to how deep surface photometry can be pushed 
without introduction of serious systematic errors. 

(ii) We find that EBL extinction effects can mimic the 
truncation of disc galaxies, but is unlikely to be the cause 
of such features in general, except maybe for the dustiest 
discs. 

(iii) EBL extinction may provide a partial explanation for 
the unexpectedly red colours reported for the halos of disc 
galaxies and hosts of local starburst galaxies. To attribute 
the colour anomalies solely to this effect would, however, re- 
quire more dust than naively expected in the faint outskirts 
of such systems. 

(iv) Based on current constraints on the dust content of 
galaxy groups and clusters, EBL extinction may lead to a 
non-negligible underestimate of intragroup and intraclus- 
ter light at the faintest surface brightness levels currently 
probed. 

(v) By combining surface photometry and direct star 
counts for systems with reasonably well-constrained opac- 
ity profiles, it should be possible to set an independent con- 
straint on the EBL. This could help resolve current contro- 
versies regarding the surface brightness of the EBL. 
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